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We compute supersymmetric contributions to pion leptonic (7172) decays in the Minimal Super- 
symmetric Standard Model (MSSM). When R-parity is conserved, the largest contributions to the 
ratio Re/fj. = F[ir + — > e + ^ e (7)]/r[7r + — > fi + ^(7)] arise from one-loop (V — A) ® (V — A) correc- 
tions. These contributions can be potentially as large as the sensitivities of upcoming experiments; if 
measured, they would imply significant bounds on the chargino and slepton sectors complementary 
to current collider limits. We also analyze R-parity violating interactions, which may produce a 
detectable deviation in i? e / M while remaining consistent with all other precision observables. 
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I. INTRODUCTION 

Low-energy precision tests provide important probes 
of new physics that arc complementary to collider 
experiments [l|, [H, In particular, effects of weak-scale 
supersymmetry (SUSY) — one of the most popular ex- 
tensions of the Standard Model (SM) — can be searched 
for in a wide variety of low-energy tests: muon (g — 2) [J , 
(3- and /j,-decay [|[ Q , parity- violating electron scatter- 
ing Q, electric dipole moment searches and SM- 
forbidden transitions like fx — > e"f Q , etc (for a recent re- 
view, see Ref. [HI). In this paper, we compute the SUSY 
contributions to pion leptonic (^12) decays and analyze 
the conditions under which they can be large enough to 
produce observable effects in the next generation of ex- 
periments. 

In particular, we consider the ratio i? e / AI , defined by 



Re/ 11 = 



T(ir + — ► e + i> e + e r v e n f) 



(1) 



The key advantage of R e u is that a variety of QCD ef- 
fects that bring large theoretical uncertainties — such as 
the pion decay constant and lepton flavor indepen- 
dent QCD radiative corrections — cancel from this ratio. 
Indeed, i? e /^ is one of few electroweak observables that 
involve hadrons and yet are precisely calculable (see [ll[ 
for discussion and Refs. [lj, EH for explicit computa- 
tions). Moreover, measurements of this quantity provide 
unique probes of deviations from lepton universality of 
the charged current (CC) weak interaction in the SM that 
are induced by loop corrections and possible extensions 
of the SM. In the present case, we focus on contribu- 
tions from SUSY that can lead to deviations from lepton 
universality. 

Currently, the two most pre cise theoretical calculations 
of R e /fj, in the SM are 



lost pre c 

0,11 



nSM 



(1.2352 ± 0.0005) x 10 
(1.2356 ± 0.0001) x 10 



(2) 



where the theoretical uncertainty comes from pion struc- 
ture effects. By utilizing chiral perturbation theory, it 
may be possible to reduce this uncertainty even fur- 
ther [l4[. Experimentally, the most precise measure- 
ments of i? e / M have been obtained at TRIUMF [l5| a nd 
PSI Taking the average of these results gives 17| 



R 



EXPT 



(1.230 ± 0.004) x 10" 



(3) 



in agreement with the SM. Future experiments at these 
facilities will make more precise measurements of R e /^, 
aimin g fo r precision at the level of < 1 x 10~ 3 (TRI- 
UMF pH) and 5 x 10~ 4 (PSI [l|). These projected un- 
certainties are close to the conservative estimate of the- 
oretical uncertainties given in Ref. (l2| . 

Deviations Ai? e / M from the SM predictions in Eq. ([2]) 
would signal the presence of new, lepton flavor-dependent 
physics. In the Minimal Supersymmetric Standard 
Model (MSSM), a non- vanishing AR^j 1 ^ may arise 
from either tree-level or one-loop corrections. In section 
HI! we consider contributions to Ai?^ sy arising from 
R-parity conserving interactions (Fig. [T]) . Although tree- 
level charged Higgs exchange can contribute to the rate 
T[ir + — > i + v(pi)\, this correction is flavor-independent 
and cancels from i? e / M . One-loop corrections induce 
both scalar and vector semileptonic dimension six four- 
fermion operators. Such interactions contribute via pseu- 
doscalar and axial vector pion decay matrix elements, 
respectively. We show that the pseudoscalar contribu- 
tions are negligible unless the ratio of the up- and down- 
type Higgs vacuum expectation values (vevs) is huge 
{v u /vd = tan/3 > 10 3 ). For smaller tan/3 the most im- 
portant effects arise from one-loop contributions to the 
axial vector amplitude, which we analyze in detail by per- 
forming a numerical scan over MSSM parameter space. 
We find that experimental observation of SUSY loop- 
induced deviations at a significant level would require 
further reductions in both the experimental error and 
theoretical SM uncertainty. Such improvements could 
lead to stringent tests of "slepton universality" of the 
charged current sector of the MSSM, for which it is often 
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FIG. 1: Representative contributions to Ai?f^ y ; (a) tree-level charged Higgs boson exchange, (b) external leg diagrams, (c) 
vertex diagrams, (d) box diagrams. Graph (a) contributes to the pseudoscalar amplitude, graphs (b,c) contribute to the axial 
vector amplitude, and graph (d) contributes to both amplitudes. 



assumed that the left-handed hrst and second gen eration 
sleptons eh and jut are degenerate (see e.g. [20| ) and 
thus AR S Y SY ~ 0. 

In section Hill we consider corrections to i? e / M from Im- 
parity violating (RPV) processes. These corrections en- 
ter at tree-level, but are suppressed by couplings whose 
strength is contrained by other measurements. In or- 
der to analyze these constraints, we perform a fit to the 
current low energy precision observables. We find that, 
at 95% C.L., the magnitude of RPV contributions to 
AR^V SY could be several times larger than the com- 
bined theoretical and anticipated experimental errors for 
the future R e /^ experiments. We summarize the main 
results and provide conclusions in section ITVl 



II. R-PARITY CONSERVING INTERACTIONS 

A. Pseudoscalar contributions 

The tree-level matrix element for ir + — > l + vg that 
arises from the (V — A) <£> (V — A) four fermion operator 



iMf v = -i2^G^V ud 



drf X PL u |tt + ) u y ^\P L vt 



2V ud F 7T G f j,mi u v P R v e , 



(4) 



where Pl.r are the left- and right-handed projection op- 
erators, 



F„ =92 A ±0.07 ±0.25 MeV 



(5) 



is the pion decay constant, is the Fermi constant ex- 
tracted from the muon lifetime, and V uc i is the (1, 1) com- 
ponent of the CKM matrix. The first error in Eq. ([5]) is 
experimental while the second arises from uncertainties 
associated with QCD effects in the one-loop SM elec- 
troweak radiative corrections to the 71^2 decay rate. The 
superscript "(0)" and subscript U AV" in Eq. (H} denote 
a tree-level, axial vector contribution. At one-loop order, 
one must subtract the radiative corrections to the muon- 
decay amplitude — since G M is obtained from the muon 
lifetime — while adding the corrections to the semilcp- 
tonic CC amplitude. The corrections to the muon-decay 



amplitude as well as lepton flavor-independent contri- 
butions to the semileptonic radiative corrections cancel 
from J? e/(tt . 

Now consider the contribution from an induced pseu- 
doscalar four fermion effective operator of the form 



ACp S = 



GpsV u d 



17(1 + 7 5 )£(i7 5 u 



(6) 



Contributions to R e /^ from operators of this form were 
considered in a model-independent operator framework 
in Ref. [2l|. In the MSSM, such an operator can arise 
at tree-level (Fig. [T^,) through charged Higgs exchange 
and at one-loop through box graphs (Fig. [TJi). These 
amplitudes determine the value oiGps- The total matrix 
element is 



iM^y+iMps = V ud F n G u m e u u (l+^)v e 



1 



Gps 



U>1 



where 



mi(m u + m d ) 



5 x 10 3 
20 



£ = e 
t = H 



(7) 



(8) 



is an enhancement factor reflecting the absence of helicity 
suppression in pseudoscalar contributions as compared to 
(V — A) g) (V — A) contributions [22j . Pseudoscalar con- 
tributions will be relevant to the interpretation of i? e / M 
if 



Gps 



G„ 



cue > 0.0005 , 



(9) 



and if Gps^e is lepton-flavor dependent. 

The tree- level pseudoscalar contribution (Fig.[T^,) gives 



G 



(0) 
PS 



mi tan (3{m u cot {3 — rrid tan 0) 



V2m 2 



(10) 



where m H + is the mass of the charged Higgs boson. 
Thus, we have 



G 



(0) 
PS 



Gn 



4 tan/3(m u cot (3 — rrid tan (3) 



(m u + md)m 2 



(11) 



H + 



3 



It is indeed possible to satisfy ^ for 

m H + 



H+ 



(12) 



Note that the combination Gpg/G^ x u>g entering Eq. ([7]) 
is independent of lepton flavor and will cancel from R e /^ . 
In principle, however, the extraction of F v from 
decay could be affected by tree-level charged Higgs ex- 
change if the correction in Eq. (J9j) is > 0.003 in magni- 
tude, corresponding to a shift comparable to the the- 
oretical SM uncertainty as estimated in Ref. [Hj]. In 
the case of charged Higgs exchange, one would require 
tan /3 > 120 [m H + /100 GeV) to generate such an effect. 

One-loop contributions to Gps are generated by box 
graphs (Fig. [TJi) . The magnitude of these contributions is 
governed by the strength of chiral symmetry breaking in 
both the quark and lepton sectors. Letting e gcncrically 
denote either a Yukawa coupling yj or a ratio m f/MsuSY 
(where / = e, [i, u, or d), we find that 



G 



(i) 

PS 



G L 



Stts^ 



( m w V 
V Msusy J 



(13) 



where the superscript "(1)" denotes one loop induced 
pseudoscalar interaction. We have verified by explicit 
computation that the 0(e) contributions vanish. The 
reason is that in each pair of incoming quarks or outgoing 
lcptons the two fermions must have opposite chirality in 
order to contribute to Gpg . Since CC interactions in the 
MSSM are purely left-handed, the chirality must change 
at least twice in each graph, with each flip generating 
a factor of e. For example, we show one pseudoscalar 
contribution in Fig. [2] that is proportional to e 2 = y^yd- 
Here, the chirality changes at the udH and v]lH vertices. 
Potentially, this particular contribution can be enhanced 
for large tan/3; however, to satisfy ([9]), we need 



tan/3 > 10 3 



(14) 



These extreme values of tan/3 can be problematic, lead- 
ing yt and y T to become nonperturbatively large. To 
avoid this scenario, we need roug hly tan < 65 (see Q 
and references therein). 

Pseudoscalar contributions can also arise through mix- 
ing of left- and right-handed scalar superpartners. Since 
each left-right mixing insertion introduces a factor of e, 

the leading contributions to Gpg will still be C(e 2 ). How- 
ever, if the triscalar SUSY-breaking parameters a/ are 
not suppressed by y/ as normally assumed, it is possible 
to have e ~ C(l), potentially leading to significant con- 
tributions. This possibility, although not experimentally 
excluded, is considered theoretically "unnatural" as it re- 
quires some fine-tuning to avoid spontaneous color and 
charge breaking (see Ref. Q for discussion). Neglecting 
this possibility and extremely large values of tan /3, we 
conclude that loop-induced pseudoscalar contributions 
are much too small to be detected at upcoming exper- 
iments. 
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FIG. 2: This contribution to Gp S is suppressed by e 2 = y^yd- 



B. Axial vector contributions 

One-loop radiative corrections also contribute to the 
axial vector matrix element. The total matrix element 
can be written as 

IMav = VudfirG^mi u v (l + j 5 )v e [1 + Arv - Ar M ] , 

(15) 

where Atv and Af^ denote one-loop contributions to the 
semileptonic and /x-decay amptlidues, respectively and 
where the hat indicates quantities renormalized in the 
modified dimensional reduction (DR) scheme. Since Af p 
cancels from R e j ill we concentrate on the SUSY contri- 
butions to Arv that do not cancel from R e /^ . It is helpful 
to distinguish various classes of contributions 



Af. 



SUSY 



A B + A GB , (16) 



where Af (Af,), 



Ay (Ay), Ab, and Aqb denote leptonic 
(hadronic) external leg (Fig.[T}3), leptonic (hadronic) ver- 
tex (Fig. [Tt), box graph (Fig. [TJi), and gauge boson 
propagator contributions, respectively. The corrections 
and Aqb cancel from R e /u, so we do not discuss 



A 9 



them further (we henceforth omit the "£" superscript). 
The explicit general formulae for Al^v,b, calculated in 
DR, are given in appendi x [Al We have verified that A^ 
and Ay agree with Ref. [23|] for case of a pure SU(2)l 
chargino/neutralino sector. 

At face value, it appears from equations (|A7HA9[) that 
carries a non-trivial dependence on MSSM pa- 

e/ ll 

rameters since the SUSY masses enter both explicitly in 
the loop functions and implicitly in the mixing matrices 
Z, defined in equations (|AHIA6[) . Nevertheless, we are 
able to identify a relatively simple dependence on the 
SUSY spectrum. 

We first consider Ai?^ sy in a limiting case obtained 
with three simplifying assumptions: (1) no flavor mixing 
among scalar superpartners; (2) no mixing between left- 
and right-handed scalar superpartners; and (3) degener- 
acy between It, and and no gaugino-Higgsino mixing. 
Our first assumption is well justified; flavor mixing in 
the slepton and squark sectors is heavily constrained by 
limits on flavor violating processes, such as /i — > e 7 0] . 

Our second assumption has minimal impact. In the 
absence of flavor mixing, the charged slepton mass matrix 



4 



decomposes into three 2x2 blocks; thus, for flavor I, the 
mass matrix in the {(-l,(-r\ basis is 



Ml 



'ir 



-±)m|cos2/3 mi (j£ — /ztan/3 
M\ — s^to| cos 2(3 



where M\ (Mj|) is the SUSY-breaking mass parameter 
for left-handed (right-handed) sleptons, ag is the coeffi- 
cient for the SUSY-breaking triscalar interaction, yg is 
the Yukawa coupling, and \i is the Higgsino mass param- 
eter. Under particular models of SUSY-breaking medi- 
ation, it is usually assumed that ag/yg ~ Msusy, and 
thus left-right mixing is negligible for the first two gener- 
ations due to the smallncss of m e and . Of course, ag 
could be significantly larger and induce significant left- 
right mixing For reasons discussed above, we neglect 
this possibility. 

We have adopted the third assumption for purely il- 
lustrative purposes; we will relax it shortly. Clearly, 
fermions of the same weak isospin doublet are not de- 
generate; their masses obey 



,2 



'ir 



Mr 



cos 2(3 + mj 
cos 2(3 + m; 



(17) 
(18) 



In addition, gaugino mixing is certainly always present, 
as the gaugino mass matrices contain off-diagonal ele- 
ments proportional to mz [see Eqs. (|A2[ IA4[) ]. However, 
the third assumption becomes valid for Msusy 3> m>z- 

Under our three assumptions, the SUSY vertex and 
external leg corrections sum to a constant that is inde- 
pendent of the superpartner masses, leading to consid- 
erable simplifications. The Bino [U(l)y gaugino] vertex 
and external leg corrections exactly cancel. The Wino 
[SU(2)i gaugino] vertex and leg corrections do not can- 
cel; rather, Ay + A^ = a/Airs^, a constant that car- 
ries no dependence on the slepton, gaugino, or Higgsino 
mass parameters. The occurrence of this constant is 
merely an artifact of our use of the DR renormaliza- 
tion scheme. (In comparison, in modified minimal sub- 
traction, we find Ay + A^ = in this same limit. 1 ') 
This dependence on renormalization scheme cancels in 
i? e //i' (In addition, this scheme-dependent constant en- 
ters into the extraction of G M ; hence, the individual decay 
widths r(-7r — > ivg) are also independent of renormaliza- 
tion scheme.) 

The reason for this simplification is that under our as- 
sumptions, we have effectively taken a limit that is equiv- 
alent to computing the one-loop corrections in the ab- 
sence of electrowcak symmetry breaking. In the limit of 
unbroken SU(2)ixU(l)y, the one-loop SUSY vertex and 



t Technically, since MS breaks SUSY, it is not the preferred renor- 
malization scheme for the MSSM. However, this aspect is not 
important in the present calculation. 



external leg corrections sum to a universal constant which 
is renormalization scheme-dependent, but renormaliza- 
tion scale-independent (23|. (For unbroken SU(2)/,, the 
SM vertex and external leg corrections yield an addi- 
tional logarithmic scale dependence; hence, the SU(2)/, 
/3-function receives contributions from both charge and 
wavefunction renormalization.) In addition, virtual Hig- 
gsino contributions are negligible, since their interac- 
tions are suppressed by small first and second generation 
Yukawa couplings. Setting all external momenta to zero 
and working in the limit of unbroken SU(2)i symmetry, 
we find that the Higgsino contributions to Al + Ay are 
2/I/327T 2 . 

In this illustrative limit, the only non-zero contribu- 
tions to AR S V SY come from two classes of box graphs 
(Fig. |TJi) — one involving purely Wino-like interactions 
and the other with both a virtual Wino and Bino. The 
sum of these graphs is 



A 



(O 



a 



12irs 2 w 



Mr 



Mf 



where we have defined 



[Fi(x L ,x Q ) + t 2 v F 2 (x B ,x L ,x Q )] 

(19) 



Fi(x Ll x Q ) 



and 



x l {xl - 2) In xl 



(xl 
xq(xq 



xq){1 - x L f 
2) Iiixq 



(20) 



{xq - x L ){l - xq) 2 (1-x l ){1-xq) 



F 2 (x b ,x l ,x q ) = - 



xb (xb + 2^/xb) In xb 



1(1 


- x B ){x B - 


- x L ){x B - x Q ) 




xl(xl + 2 y 


/xb ) In xl 


(1- 


- x l ){xl - 


xb)(xl - xq) 




xq(xq + 2 


s/xb) In Xq 


(1- 


-x Q )(x Q - 


xl){x q - X B )_ 



(21) 



where xb = Mf/Mf , xl = m~/M|, and xq 



-JMl 



and m,Q of the Bino, Wino, left- 



with masses Mi, M 2 , mj 
handed ^-flavored slepton, and left-handed 1st generation 
squark, respectively. Numerically, we find that always 
Fi ^s> F 2 ; the reason is that the sum of Bino- Wino graphs 
tend to cancel, while the sum of pure Wino graphs all 
add coherently. Hence, Bino exchange (through which 
the term proportional to F 2 arises) does not significantly 
contribute to ARfV SY . 

In Fig. [31 we show Fi(xl, xq) as a function of xl for 
fixed xq. Since F\ is symmetric under xl <-» xq, Fig. [3J 
also shows F\ as a function of xq , and hence how A b de- 
pends on m^ L . For xl, xq ~ 1, we have F\ ~ 0(1), while 
if either ij, > 1 or > 1, then F\ — > 0, which corre- 
sponds to the decoupling of heavy sleptons or squarks. 
There is no enhancement of for xl <C 1 or xq <C 1 
(i.e. if M 2 is very heavy) due to the overall 1/M| sup- 
pression in (TT9"]) . 
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FIG. 3: The box graph loop function F\(xl,xq) as a function 
of xl = m~/M;f for several values of xq = m^/M?. For 
xl ~ xq ~ 1 (i.e. SUSY masses degenerate), Fi(xl,xq) ~ 
1. For xl > 1 or xq 3> 1 (i.e. very massive sleptons or 
squarks), Fi(xl,xq) — > 0. 



FIG. 4: ARf^ Y versus fi, with fixed parameters Mi = 100 
GeV, M 2 = 150 GeV, m ?L = 100 GeV, m fe = 500 GeV, 
m; t = 200 GeV. Thin solid line denotes contributions from 
(Ay + Al) only; dashed line denotes contributions from Ab 
only; thick solid line shows the sum of both contributions to 



The total box graph contribution is 



AR s y sY 



2 Rc[A^ - a£ } ] 



( m w X 
\ M 2 J 



6tts^ 



(22) 



Mf M| 



Mr 



M| 



Clearly Ai?^ SY vanishes if both sleptons are degenerate 
and is largest when they are far from degeneracy, such 
that mj L 3> nnji L or mg L <C Trvp L . In the latter case, we 
have 



Ai? SUSY 



pSM 

e / A 1 



< 0.001 x 



100 GeV 

Msusy 



(23) 



for e.g. Msusy = M 2 ~ m Si ~ itIq l < mjj L . 

We now relax our third assumption to allow for 
gaugino-Higgsino mixing and non-degeneracy of I and 
Vi. Both of these effects tend to spoil the universality of 
Ay + A L , giving 



4tts^ 



/- 0.001/ 



(24) 



The factor / measures the departure of Ay + A^ from 
universality. If the SUSY spectrum is such that our 
third assumption is valid, we expect / — > . For re- 
alistic values of the SUSY parameters, two effects lead 
to a non-vanishing /: (a) splitting between the masses 
of the charged and neutral left-handed sleptons that re- 
sults from breaking of SU(2)l, and (b) gaugino-Higgsino 
mixing. The former effect is typically negligible. To see 



why, we recall from Eq. (fTT|) that 



1 + 



(25) 



where we have neglected the small non-degeneracy pro- 
portional to the square of the lepton Yukawa coupling. 
We find that the leading contribution to / from this non- 
degeneracy is at least 0(m^/m~), which is < 0.1 for 
mj> 2M W . 

Significant gaugino mixing can induce / ~ 0(1)- The 
crucial point is that the size of / from gaugino mixing is 
governed by the size of M 2 . If M 2 ^> mz, then the Wino 
decouples from the Bino and Higgsino, and contributions 
to Ay + Al approach the case of unbroken SU(2)l. On 
the other hand, if M 2 ~ mz, then Ay + can differ 
substantially from a/Aits^- 

In the limit that mj 3> M2 ( 
a decoupling scenario where 
and thus / = 0. Hence, a significant contribution to 
Ai? e / M requires at least one light slepton. However, re- 
gardless of the magnitude of /, if 
corrections will cancel from R. 



= e, pi), we also have 



m^ L , then these 



//'■ 



It is instructive to consider the dependence of individ- 
ual contributions and Ay+A^ to AR~¥f , as shown 
in Figs. 2] and [5j In Fig. [3J we plot the various contribu- 
tions as a function of /x, with Mi = 100 GeV, M 2 = 150 
GeV, TOe L = 100 GeV, m ?i = 500 GeV, m 5i = 200 GeV. 
We see that the Ay + A^ contributions (thin solid line) 
vanish for large fx, since in this regime gaugino-Higgsino 
mixing is suppressed and there is no Ay + A^ contribu- 
tion to AR^Vf Y . However, the A# contribution (dashed 
line) is nearly /^-independent, since box graphs with Hig- 
gsino exchange (which depend on [i) are suppressed in 
comparison to those with only gaugino exchange. In Fig. 



G 




FIG. 5: AR^ Y /R^ as a function of M 2 , with /j, = 200 
GeV and all other parameters fixed as in Fig. [7J Each line 
shows the contribution indicated as in the caption of Fig. fj| 



[5j we plot these contributions as a function of M2, with 
/j = 200 GeV and all other parameters fixed as above. 
We see that both Ay + At, and A b contributions vanish 
for large M2. 

One general feature observed from these plots is that 
Ay + Al and As contributions tend to cancel one 
another; therefore, the largest total contribution to 
ARfVf Y occurs when either Ay + Al or A^ is sup- 
pressed in comparison to the other. This can occur in the 
following ways: (1) if /1 3> mz, then Ab may be large, 
while Ay + A^ is suppressed, and (2) if mz L , 3> mz, 
then Ay + Al may be large, while A# is suppressed. In 
Fig. O we have chosen parameters for which there is a 
large cancellation between Ay + Al and Ab- However, 
by taking the limits /1 — ► 00 or mz L , m ^ L ~ * ^-^eAj 
would coincide the As or Ay + A^ contributions, respec- 
tively. 

Because the Ay + Al and A^ contributions tend to 
cancel, it is impossible to determine whether or JIl is 
heavier from R e /^ measurements alone. For example, a 
positive deviation in R e / I1 can result from two scenarios: 
(1) ARfj 1 ^ is dominated by box graph contributions 
with rrie L < mp L , or (2) AR^ SY is dominated by Ay + 
A l contributions with mz L > m^ L . 

Guided by the preceding analysis, we expect for 



• The maximum contribution is 
0.001. 



AR^/R 



< //' 



• Both the vertex + leg and box contributions are 
largest if M2 ~ 0(mz) and vanish if M2 3> mz- If 
M2 ~ 0(mz), then at least one chargino must be 
light. 

• The contributions to AR S ¥ SY vanish if m? r = 



and are largest if either m 



Mi < m ?L 



0.0002 



0.1 



10 



m g/ m^ 



mp L » rrie L . 



FIG. 6: AR^^ Y as a function of the ratio m^ L jm^ L . Pa- 
rameter points which obey the LEP II bound are dark blue; pa- 
rameter points which violate the LEP II bound are light blue. 



• The contributions to Ai?^ sy are largest if ei or 
J1 L is 0(m z ). 

• If /j > mz, then the lack of gaugino-Higgsino 
mixing suppresses the Ay + Al contributions to 

• If m% L , m^ S> mz, then the contributions to 

are suppressed due to squark decoupling. 

• If ul, cIl, and /i are all 0{mz), then there may be 
cancellations between the Ay + Al and Ab con- 
tributions. ARfVf Y is largest if it is dominated by 
either Ay + Al or Ab contributions. 

We now study AR^^ Y quantitatively by making a 
numerical scan over MSSM parameter space, using the 
following ranges: 

mz/2 < {Mi, \M 2 \,\n\, mz L } < 1 TcV 

m z /2 < {m^,mzj < 5 TeV (26) 
1 < tan/3 < 50 
sign(/x), sign(M 2 ) = ±1 , 

where m^ L , mji L , and m t are determined from Eqs. 

(EES). 

Direct collider searches impose some constraints on 
the parameter space. Although the detailed nature of 
these constraints depend on the adoption of various as- 
sumptions and on interdependences on the nature of the 
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FIG. 7: AT?; 



/M 



a function of Min[m-g L , mjx L ], the 



of the lightest first or second generation charged slepton. Pa- 
rameter points which obey the LEP II bound are dark blue; pa- 
rameter points which violate the LEP II bound are light blue. 



FIG. 8: 



versus m xi , the mass of the lightest 



chargino. Parameter points which obey the LEP II bound are 
dark blue; parameter points which violate the LEP II bound 
are light blue. 



MSSM and its spectrum [lTJ, we implement them in a 
coarse way in order to identify the general trends in cor- 
rections to Re/fi- First, we include only parameter points 
in which there are no SUSY masses lighter than mz/2. 
(However, the current bound on the mass of lightest neu- 
tralino is even weaker than this.) Second, parameter 
points which have no charged SUSY particles lighter than 
103 GeV are said to satisfy the "LEP II bound." (This 
bound may also be weaker in particular regions of pa- 
rameter space.) 

Additional constraints arise from precision electroweak 
data. We consider only MSSM parameter points whose 
contributions to oblique parameters S, T, and U agree 
with Z-pole measurements at 95% CL [131 ■ Because 
we have neglected the 3rd generation and right-handed 
scalar sectors in our analysis and parameter scan, we do 
not calculate the entire MSSM contributions to S, T, 
and U. Rather, we only include those from charginos, 
neutralinos, and 1st generation left-handed scalar super- 
partners. Although incomplete, this serves as a conser- 
vative lower bound; in general, the contributions to S, 
T, and U from the remaining scalar superpartners (that 
we neglect) only causes further deviations from the mea- 
sured values of the oblique parameters. In addition, we 
include only the lightest CP-even Higgs boson with mass 
rrih = 114.4 GeV, neglecting the typically small contri- 
butions from the remaining heavier Higgs bosons. 

We do not impose other electroweak constraints in 
the present study, but note that they will generally lead 
to further restrictions. For example, the results of the 



E821 measurement of the muon anomalous magnetic mo- 
ment [24| tend to favor a positive sign for the /i param- 
eter and relatively large values of tan (3. Eliminating the 
points with sign(/z) = — 1 will exclude half the parameter 
space in our scan, but the general trends are unaffected. 

We show the results of our numerical scan in Figs. [5]- 
[5] Parameter points which satisfy the LEP II bound are 
dark blue; those which do not arc light blue. In Fig. [HI wc 
show ARfj^f Y /Rg/fj, as a function of the ratio of slepton 

masses me L /mp L . If both slcptons are degenerate, then 
AR S V SY vanishes. Assuming the LEP II bound, in order 

e/ fj. ° 3 

for a deviation in i? e /^ to match the target precision at 
upcoming experiments, we must have 



5R, 



e/fj. 



AR™ SY /R ( 



//' 



> 0.0005 



(27) 



and thus m^ L /m-p L > 2 or mi 



■Jms L >2. 



In Fig. wc show ARfZ SY /R as a function of 

M.m[mz L , mpL,], the mass lightest first or second gen- 
eration slepton. If the lighter slepton is extremely heavy, 
then both heavy sleptons decouple, causing AR^Vf Y to 
vanish. Assuming the LEP II bound, to satisfy (|2T|) . we 
must have rrie L < 300 GeV or m^ L < 300 GeV. 

In Fig. [SJ we show AR^Vf Y /R,,/^ as a function of m x i, 

the lightest chargino mass. If m x i is large, AR^¥f Y van- 
ishes because M% must be large as well, suppressing Ab 
and forcing Ay and Al to sum to the flavor independent 
constant discussed above. Assuming the LEP II bound, 
to satisfy (|2"T)) . we must have m x i < 250 GeV. 
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FIG. 10: Tree-level RPV contributions to R, 



governed by both the sfermion masses and by the pa 
rameters A' llfe 
interactions: 



and W 21k that are the coefficients in RPV 



&RPV, AL=1 — Kjk^iQj^k 

Defining 0, ^ 



AU(/) = 



A; 



ijk I 



4V2G,, 



>0, 



(28) 



(29) 



FIG 



9: Contours indicate the largest values of 8R^j'^ Y ob- 



tained by our numerical parameter scan 126)1 . as a function 
of and m; t . The solid shaded regions correspond to the 
largest values of 5R.f¥^ Y within the ranges indicated. All val- 
ues of 5ReJ^ Y correspond to parameter points which satisfy 
the LEP II bound. 



Finally, we illustrate the interplay between Ay + Al 
and Ab by considering SR^Vf Y as a function of and 

m„ L ■ In Fig.[9l we show the largest values of SR^Vf Y ob- 
tained in our numerical parameter scan, restricting to pa- 
rameter points which satisfy the LEP II bound. The solid 
shaded areas correspond to regions of the |/i|-ms L plane 
where the largest value of 5R S ¥ SY lies within the indi- 

° e/n 

catcd ranges. It is clear that 5RfVf Y can be largest in 
the regions where either (1) \x is small, m,u L is large, and 



the largest contributions to AR 



SUSY 



are from Ay + At, , 



or (2) /i is large, ra^ L is small, and the largest contri- 
bution to ARSys Y is from Ab- If both /i and m„ £ arc 

light, then Ai?^ t sy can still be very small due to can- 
cellations, even though both Ay + A^ and Ab contri- 
butions are large individually. More precisely, to satisfy 
p7|) . we need either fi < 250 GeV, or fi > 300 GcV and 



rriu L < 200 GeV. 



III. CONTRIBUTIONS FROM R- PARITY 
VIOLATING PROCESSES 



In the presence of RPV interactions, tree-level ex- 
changes of sfermions (shown in Fig. IT0|) . lead to viola- 
tions of lepton universality and non-vanishing effects in 
R e /^ ■ The magnitude of these tree-level contributions is 



contributions to i? e / M from RPV interactions are 
AR R , PV 

- OA' _ OA' 
R SM - ZZA llfe Zl ±21k- 

e/fj, 



(30) 



Note that RPV contribution to the muon lifetime (and, 
thus, the Fermi constant G^) cancels in i? e / M , therefore 
does not enter the expression. 

The quantities A'- k etc. are constrained by existing 
precision measurements and rare decays. A summary 
of the low energy constraints is given in Table III of 
Ref. , which includes tests of CKM unitarity (primar- 
ily through RPV effects in superallowed nuclear /3-decay 
that yields a precise value of \V u d\ [HI), atomic parity 
viola ting (PV) measurements of the cesium weak charge 
Qw EM j the ratio i? e //i itself [H, 0j[ , a comparison of 
the Fermi constant G M with the appropriate combination 
of a, mz, and sin 2 6\y [29| . and the electron weak charge 
determined from SLAC E158 measurement of parity vi- 
olating M0ller scattering [3Cj. 

In Fig. [TT]we show the present 95% C.L. constraints 
on the quantities A' llk and A' 21k obtained from the afore- 
mentioned observables (interior of the blue curve). Since 
the A'- k are positive semidefinite quantities, only the re- 
gion of the contour in the upper right hand quadrant are 
shown. The green curve indicates the possible implica- 
tion of a future measurement of the proton weak charge 
planned at Jefferson Lab [3l|, assuming agreement with 
the Standard Model prediction for this quantity and the 
anticipated experimental uncertainty. The dashed red 
curve shows the possible impact of future measurements 
of Rf./fi, assuming agreement with the present central 
value but an overall error reduced to the level antici- 
pated in Ref. [H]; with the error anticipated in Ref. [l9| 
the width of the band would be a factor of two smaller 
than shown. 
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FIG. 11: Present 95% C.L. constraints on RPV parameters 
A'jiki 3 ' = 1)2 that enter i? e / M obtained from a fit to precision 
electroweak observables. Interior of the dark blue contour 
corresponds to the fit using the current value of AJ? e/ / JJ /7?f/^f 



[151 . 116j , while the dashed red contour corresponds to the fit 
using the future expected experimental precision [3, assum- 
ing the same central value. The light green curve indicates 
prospective impact of a future measurement of the proton 
weak charge at Jefferson Lab [3l| . 



IV. CONCLUSIONS 

Given the prospect of two new studies of lepton uni- 
versality in tt£ 2 decays 0, with experimental errors 
that are substantially smaller than for existing measure- 
ments and possibly approaching the 5 x 10~ 4 level, an 
analysis of the possible implications for supersymmetry 
is a timely exercise. In this study, we have considered 
SUSY effects on the ratio R e/fl in the MSSM both with 
and without R-parity violation. Our results indicate that 
in the R-parity conserving case, effects from SUSY loops 
can be of order the planned experimental error in par- 
ticular, limited regions of the MSSM parameter space. 
Specifically, we find that a deviation in i? e / AI due to the 
MSSM at the level of 



0.0005 < 



ar susy 



R, 



/!■> 



< 0.001 , 



(31) 



implies (1) the lightest chargino %i is sufficiently light 

m xl < 250 GeV , 

(2) the left-handed selectron ei and smuon ju^ are highly 
non-degenerate : 



> 2 



< 



1 



(3) at least one of or /i£ must be light, such that 

< 300 GeV, 



m e L ~ 300 GeV or m? 



Two general observations emerge from Fig. [TTJ First, 



values of A' 21k 



and A 



life 



given the present constraints, 
differing substantially from zero are allowed. For val- 
ues of these quantities inside the blue contour, ARfV, SY 
could differ from zero by up to five standard deviations 
for the error anticipated in Ref. Such RPV effects 

could, thus, be considerably larger than the SUSY loop 
corrections discussed above. On the other hand, agree- 
ment of Re/fj, with the SM would lead to considerable 
tightening of the constraints on this scenario, particu- 
larly in the case of A' 21fc , which is currently constrained 
only by i? e / M and deep inelastic v {y) scattering (32J. 



The presence of RPV interactions would have signifi- 
cant implications for both neutrino physics and cosmol- 
ogy. It has long been known, for example, that the exis- 
tence of AL = ±1 interactions — such as those that could 



enter R, 



will induce a Majorana neutrino mass [3c 



while the presence of non- vanishing RPV couplings would 
imply that the lightest supcrsymmetric particle is unsta- 
ble and, therefore, not a viable candidate for cold dark 
matter. The future measurements of R e /^ could lead 
to substantially tighter constraints on these possibilities 
or uncover a possible indication of RPV effects. In ad- 
dition, we note that the present uncertainty associated 
with RPV effects entering the ir^ decay rate would af- 
fect the value of F n at a level about half the theoretical 
SM uncertainty as estimated by Ref. [l2j ■ 



and (4) the Higgsino mass parameter fi and left-handed 
up squark mass m« L satisfy either 

< 250 GeV 



or 



|m| > 300 GeV, m UL < 200 GeV. 



Under these conditions, the magnitude Ai?^ sy may fall 
within the sensitivity of the new R e /^ measurements. 

In conventional scenarios for SUSY-brcaking media- 
tion, one expects the left-handed slepton masses to be 
comparable, implying no substantial corrections to SM 
predictions for i? e / M . Significant reductions in both ex- 
perimental error and theoretical, hadronic physics uncer- 
tainties in R s f [ would be needed to make this ratio an 
effective probe of the supcrpartner spectrum. 

On the other hand, constraints from existing precision 
electroweak measurements leave considerable latitude for 
observable effects from tree-level superpartner exchange 
in the presence of RPV interactions. The existence of 
such effects would have important consequences for both 
neutrino physics and cosmology, as the presence of the 
AL RPV interactions would induce a Majorana mass 
term for the neutrino and allow the lightest supcrpartner 
to decay to SM particles too rapidly to make it a viable 
dark matter candidate. Agreement between the results 



of the new R 



e/ii 



measurements with RfJ 1 ^ could yield 



significant new constraints on these possibilities. 
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APPENDIX A: GENERAL RADIATIVE CORRECTIONS IN THE MSSM 

The MSSM Lagrangian and Feynman rules [34j are expressed in terms of chargino and neutralino mixing matrices 
Z± and Zm, respectively, which diagonalize the supcrpartner mass matrices, defined as follows. The four neutralino 
mass eigenstates Xi are related to the gauge eigenstates ip° = (B, W 3 ,H®, H®) by 



Vi — X-3 > 



where 



7 T 



/ Mi 
M 2 
-Cp s w m z Cp c w m z 
V sp swmz —sp cwmz 



-cp s w m z 
cp cwm z 





J N 



( m xl 




V o 











to^o 




(Al) 



(A2) 



l) 



is the diagonalized neutralino mass matrix. The chargino mass eigenstates xt are related to the gauge eigenstates 
i>+ = (W+,H+) and V" = (W~,H-) by 



7 ZJ -v d 



where 



z t ( M 2 
\ \[2cpm w 



\plspmw 



"Xl 





(A3) 



(A4) 



is the diagonalized chargino mass matrix. We note that the off-diagonal elements which mix gauginos and Higgsinos 
stem solely from electroweak symmetry breaking. 

The charged slcpton mass eigenstates Li are related to the gauge eigenstates £ = ifiL^L^h^R^R^TB) by 



- 7 % 3 



Z\ J L 



3 1 



where 



2 

TO- 



Z\ Mf Z L = 



(A5) 



(A6) 







is the diagonalized slepton mass matrix. There are two classes of off-diagonal elements in M~ which can contribute to 
slcpton mixing: mixing between flavors and mixing between left- and right-handed components of a given flavor, both 
of which arise through SUSY-breaking terms. (Left-right mixing due to SUSY-preserving terms will be suppressed by 
Tni/mj and is irrelevant for the first two generations.) 

Similarly, up- type squarks, down- type squarks, and sneutrinos have mixing matrices Zjj, Zd, and Z v , respectively, 
defined identically to Zl — except for the fact that there are no right-handed sneutrinos in the MSSM and thus there 
are only three sneutrino mass eigenstates. 

There are three types of contributions to Ai?^ sy in the MSSM: external leg, vertex, and box graph radiative 
corrections. The leptonic external leg corrections (Fig. [TJa) are 



16tts w 



\Z$t w - Z 2 ^\ 2 B(m x o,m^) + 2 |Zi fe | 2 B{m Xk ,m z ) 
\Z%t w + Z%\ 2 B{m x o,m Zi ) + 2 \Z\ k \ 2 B{m Xk ,m^) 



(A7) 
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where the loop function is [H| 

B(mi, m 2 ) = / dx x In — 5 

Jo \ m i 

The leptonic vertex corrections (Fig. []};) are 



A/ 2 



(1 — x) + m\x 



A 



(/) 



(A8) 



(z#* - t w z|* f (z^z| fe - ^^+ fe ) c 



(Z^*Z U ' + -^Z^*Z 2fc ) m x om Xfc Ci(m x o,r 



1 



2 (Z# + % Z#) Zl k [ (Z 2 J*Z}_ k + -j=Z 3 J*Z 2 _ k ) C 2 (m Xk ,m Zt ,m x o) 



V2 



Z^Zf) m x om Xk C 1 (m Xh ,m z ,,m x a) J , 



with loop functions 



Ci (mi, 7712,1713) = / dxdy — 5- 

7o 



+ m|y + 7713(1 — x — y) 



C2 (mi, 7712, 7713) = / cfa; In 
Jo 

The corrections from box graphs (Fig. [TJi) are 



M 2 

m 2 x + + 7773(1 — x — y) 



A 



U) 



Mr 



( |Z u f (Z^"» + tw^ m *) {Z 2 N m - l -t w Z]l n ) D x {m xam ,m Ih ,m Xh ,m Zi ) 



(A9) 



+ \Z^\ 2 (Z 2 N m ~ t W Zj^) (Z% m * + -twZ 1 ^*) Di (777^,777^, 777^, 777^) 

+ Zi^' (Z 2 ,™ - t w Z#") (Z 2 ?" - i%Z#") m x o m m X: ,. £2(777^, m^, m Xj ,m Vi ) 

+ Z lk Z lk (Z 2 ,™* +t w Z 1 N m *) (Zjj n * + \t w Z 1 ^*) m x om Xk D^m^m^m^m^) ) , 



with loop functions 



D„(mi, 777 2 , 7773, m 4 ) 



dx dy dz . „ ^ ^ o ? * m 
[mfa; + 7772?; + + m|(l — x — y — z)J™ 



In formulae (|A7IIA9|) , 7=1 corresponds to ir — > e z/ e and 1 = 2 corresponds to 7r — > /j All other indeces 
are summed over. We use DR renormalization at scale M. We have defined ivr = tan#jy and sw = singly. We 
have neglected terms proportional to either Yukawa couplings or external momenta (which will be suppressed by 
0(m v /MsuSY))- Finally, the SUSY contribution to i? e / M is 



Afl SUSY 



R, 



= 2 Re[A« - A<? + A« - A« + A« - A^] 



(A10) 



In addition, the following are some useful formulae needed to show the cancellations of vertex and leg corrections 
in the limit of no supcrpartner mixing: 

62(7711, 777 2 , 777i ) = B(m 2 ,mi) 
2m\ Ci(mi, 7772 , 777i) — 2 B(mi, m 2 ) + 2 B(m 2 , mi) = 1 . 
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